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Abstract 

Identification of tlie unknown parameters and orders of fractional chaotic sys- 
tems is of vital significance in controlling and synchronization of fractional- 
order chaotic systems. However there exist basic hypotheses in traditional 
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estimation methods, that is, the parameters and fractional orders are par- 
tially known or the known data series coincide with definite forms of frac- 
tional chaotic differential equations except some uncertain parameters and 
fractional orders. What should I do when these hypotheses do not exist? 

In this paper, a non-Lyapunov novel approach with a novel united mathe- 
matical model is proposed to reconstruct fractional chaotic systems, through 
the fractional-order differential equations self-growing mechanism by some 
genetic operations ideas independent of these hypotheses. And the cases 
of identifying the unknown parameters and fractional orders of fractional 
chaotic systems can be thought as special cases of the proposed united 
mathematical reconstruction method in non-Lyapunov way. The problems 
of fractional-order chaos reconstruction are converted into a multiple modal 
non-negative special functions' minimization through a proper translation, 
which takes fractional-order differential equations as its particular indepen- 
dent variables instead of the unknown parameters and fractional orders. And 
the objective is to find best form of fractional-order differential equations such 
that the objective function is minimized. Simulations are done to reconstruct 
a series of hyper and normal fractional chaotic systems. The experiments' 
results show that the proposed self-growing mechanism of fractional-order 
differential equations with genetic operations is a successful methods for 
fractional-order chaotic systems' reconstruction, with the advantages of high 
precision and robustness. 

Keywords: Fractional-order chaotic system, reconstruction, 
non-Lyapunov, non-negative special functions' minimization, fractional 
differential equations as independent variables 
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1. Introduction 

Although the concept of fractional calculus was proposed by Leibniz 300 
years ago [1-5], it become a powerful tool only during the past decades to de- 
scribe the dynamics of complex systems in area of science and engineering [6- 
11]. 

The applications of fractional differential equations began to appeal to 
related scientists[12-37] after the symbolic event that Mandelbort discovered 
that there were a lot of fractional dimension phenomena in nature [38] in fol- 
lowing areas, bifurcation, hyperchaos, proper and improper fractional-order 
chaos systems and chaos synchronization [6-37]. And the integer order chaos 
systems are the special cases of fractional order chaos systems. It is much 
more important to research the hyper fractional-order chaos systems, both 
proper and improper fractional orders. 

For the fractional-order chaos systems, most of the control and syn- 
chronization methods are generalized conclusive but more complicated ap- 
proaches from the methods for normal and hyper chaos system. And there 
are no results on hyper or improper fractional-order chaos systems up to now. 

However, there are some systematic parameters and orders are unknown 
for the fractional-order chaos systems in controlling and synchronization. It 
is difficult to identify the parameters in the fractional-order chaotic systems 
with unknown parameters. 

Hitherto, there have been at two main approaches in parameters' identi- 
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fication for fractional-order chaos systems. 

The first is a Lyapunov way using synchronisation. And there have been 
few resuhs on parameter estimation method of fractional-order chaotic sys- 
tems based on chaos synchronization [39]. However, the design of controller 
and the updating law of parameter identification is a task with technique 
and sensitively depends on the considered systems. 

The second is a non-Lyapunov way using some artificial intelligence meth- 
ods, such as differential evolution[17] and particle swarm optimization[19]. 
Although there are many methods in parameters estimation for integer order 
chaos systems, to the best of our knowledge, a little work in non-Lyapunov 
way has been done to the parameters and orders estimation of fractional- 
order chaos systems [6, 17, 19]. 

For example, we consider the following fractional-order chaos system. 

^D',Y{t) = f{Y{t),Yo{t),e) (1) 

where Y(t) = {yi{t),y2{t), ...,yn(t))'^ G 3?" denotes the state vector. 9 = 
{9i, 02, 9n)^ denotes the original parameters, q = (gi, q2, qn), (0 < gj < 
1,2 = 1, 2, n) is the fractional derivative orders. 

f{Yit),Yo{t),e)= if,j2,...Jn)\^Yit),Yoim 

Normally the function / is known. And the 6, q are unknown , then the 
6 = (6'i, 02, 9n, qi, q2, qn) will be the parameters to be estimated. 
Then a correspondent system are constructed as following. 

^D',Y{t) = f{Y{t),Yo{t)J) (2) 
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where Y{t), 9, q are the correspondent variables to those in equation (1), and 
function / are the same. The two systems (1) (2) have the same initial 
condition Yo(t). 

And the objective function to be minimized by artificial intelligence meth- 
ods can be defined as following. 

N 



i=l 



G 



1 ^ II 



i=l 
N 



1=1 



1 ^ II 

i=l 



Then the objective is obtained as following, 



9* = argminF 



(3) 



(4) 



(5) 



(6) 



(7) 



where F in equation (7) can be any one of equations (3), (4), (5), (6). 

However, there exist basic hypotheses in traditional non-Lyapunov esti- 
mation methods [6, 17, 19]. That is, the parameters and fractional orders 
are partially known or the known data series coincide with definite forms 
/ = (/i) •••5 /n) of fractional chaotic differential equations except some 
uncertain parameters and fractional orders 6 = {9i,92, ■■.,9n,qi,q2, ■■■,qn)■ 
Wh.eit should 1 do when these hypotheses do not exist? That is, when some 
the fractional chaotic differential equations / = (/i, /2, /„) are unknown, 
how to identify the fractional system? That is. 



(/i,/2,...,/n)* = arg min F (8) 

(/l,/2v,/n) 

Now the problem of parameters estimation (7) become another much 
more comphcated question, to find the forms of fractional order equations as 
in (8). In another word, it is fractional-order chaos reconstruction problem 
now. 

However, to the best of authors' knowledge, there are no methods in non- 
Lyapunov way for fractional order chaotic systems' reconstruction so far. 
The objective of this work is to present a novel simple but effective approach 
to reconstruct hyper, proper and improper fractional chaotic systems in a 
non-Lyapunov way. In which, an idea of self growing equations through 
only genetic operations are used to identify unknown forms of differential 
equations / = (/i, /2, /„), in which no extra evolutionary algorithms to 
be added or be compounded. And the illustrative reconstruction simulations 
in different chaos systems system are discussed respectively. 

The rest is organized as follows. Section 2 give a simple review on non- 
Lyapunove parameters estimation methods for fractional-order chaos sys- 
tems. In Section 3, a novel united mathematical model for fractional chaos 
reconstruction are proposed in a non-Lyapunov way. Then all the existing 
non-Lyapunov methods in Section 2 become special cases of the new united 
model, such as non-Lyapunov parameters identification for normal and frac- 
tional chaos systems, non-lyapunov reconstruction methods for normal chaos 
system. In section 4 a novel methods with proposed united model in Section 
3 is introduced to reconstruct the hyper, improper and normal fractional 
chaos systems. And simulations are done to a series of different fractional 
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chaos systems and their special cases in three dimension. Conclusions are 
summarized briefly in Section 5. 

2. Non-Lyapunove reconstruction methods for normal chaos sys- 
tems 

To solve the question proposed in Section 1, a new ideas are introduced to 
generalize the normal parameters' estimation problems into fractional order 
chaos reconstruction problems. In the sense of the differential equations for 
chaos system reconstruction, the parameters' identification for normal chaos 
system cases[40-51] and the for the fractional Lorenz system[13, 18, 34] [26, 
52] can be thought as special cases of fractional order chaos reconstruction, 
when the exact forms of fractional order chaotic differential equations / = 
(/i, /2, /„) are available but some parameters unknown. 

To propose a novel reconstruction way for fractional-order chaos sys- 
tem, we briefly reviewed the normal chaos reconstruction methods as follow- 
ing. Firstly, we introduce the basic artiflcial intelligent symbolic regression 
method Genetic programming (GP) . Secondly, the non-Lyapunov way for 
normal chaos reconstruction. 

2.0.1. Genetic Programming 

To have a deep understanding of reconstruction method in non-Lyapunov 
way, flrstly we introduce some basic ideas of GP. 

GP proposed by J. R. Koza in the early 1990s[53, 54], belongs to evo- 
lutionary algorithms, used to automatically generate the complex function 
between the input and output data. The basic idea of GP is similar to 
genetic algorithm. However the more complexity of structure developed in 
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GP's self-adaptive "Program" evolution simulation in finding the best forms 
of the functions. Because GP uses a much more natural and fiexible repre- 
sentation, then it has a wide applications in fields of artificial intelligence, 
machine learning, control, molecular biology and etc. 

According to Darwinian principle of natural selection, individuals with 
high fitness survived by genetic operators, including crossover, mutation, 
gene duplication, gene deletion, to the next generation of groups. As the 
evolution continues, the best or the satisfied optimal solution of the given 
problem will appears. There are five major preparatory steps for the basic 
GP to be specified to be evolved. 

• Terminals set. The independent variables, zero- argument functions and 
random constants. 

• Set of primitive functions. For instances, +, — , x, ^, sin, cos, log, ... 
combine in term of tree structure or the other structure for the in- 
dividuals. 

• Fitness measure. (Explicit or implicit methods measuring the fitness 
of individuals. Similar to the genetic algorithm, an evaluation criteria 
is then defined as individuals's fitness values. 

• Controlling parameters for the GP's evolution. 

• Termination criterion and method for designating the result. 

2.0.2. Evolutionary Reconstruction methods 

The evolutionary reconstruction methods here can be concluded into two 
categories. 
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The first is a compound symbolic regression(SR) method[55] with GA, 
differential evolutionary algorithms or self-organizing migrating algorithm 
(SOMA) [56-59] and the other evolutionary algorithms. However it can be 
used for reconstruction of only one differential equation of chaotic systems 
such as for Lorenz and the other systems[58, 60]. The differences between 
Ref. [59] and the series researches of Zelinka[56-58] are the fix candidate 
coding table pre-defined [59] (that is, the basic components of the function 
form to be searched) and fitness evaluation through Lyapunov exponents. 

The second is with model with the parameters unknown in a co-evolution. 
And the models are to be reconstructed by GP. The optimal parameters' 
combinations of the reconstructed model are identified by GA, PSO[61-66]. 

It is also noticed that there is no method for reconstruction for hyperchaos 
systems and fractional chaos systems so far. 

3. Unified Novel Mathematical model for fractional chaos recon- 
struction in non-Lyapunov way 

Unlike the normal Lyapunov methods [67-69], for the normal chaos sys- 
tems, there have been some results on chaos reconstruction for / = (/i, /2, fn) 
with the non-Lyapunove methods mainly from symbolic regression through 
genetic programming(GP)[53, 54, 70], and some evolutionary algorithms[56, 
58, 60, 62-66, 71-73]. 

However, to the best of authors' knowledge, there are no direct methods 
for fractional chaotic systems' reconstruction. 

we take the fractional Lorenz system[13, 18, 34] for instance, which is 
generalized from the first canonical chaotic attractor found in 1963, Lorenz 
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system [74]. 

a,Bfx = a ■ - x); 

(9) 

aD^z = X ■ y — b ■ z. 
L = {x,y,z) 

where gi, g2, Q'3 are the fractional orders. When (gi, g2, Q's) = (0.993,0.993,0.993), 
a = 10,7 = 28,6 = 8/3, a = 0, intimal point (0.1,0.1,0.1) system (9) is 
chaotic. Generally when the dimension 

5^ = gi + 52 + gs > 2.91 

for fractional system (9) is chaotic. 

The original system is Eq. (9) or matrix data observed from the actual 
system. Then we generate the three unknown equations fi{x, y, z),i = 1,2,3 
for reconstructing the chaos system (10). 

aD^'x = fi{x, y,~z) = a ■ {y - x); 
aDTy = f2{x,y,z)=-i-x-x-z-y; 

(10) 

aDl^'z = fs{x,y,z) = x-y-b-z. 
L = {x,y,z) 

Let f-s{x,y,z) be unknown in system (10), and the objective here is to 
reconstruct the /a. 

And there exist several definitions of fractional derivatives. Among these, 
the Griinwald-Letnikov (G-L), Riemann-Liouville (R-L) and the Caputo frac- 
tional derivatives are the commonly used[2-5, 75, 76]. And G-L, R-L and 
Caputo fractional derivatives are equivalent under some conditions [77]. 
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The continuous integro-differential operator[78, 79] is defined as 



l,g = 0; 

We consider the continuous function f{t). The G-L fractional derivatives 
are defined as following. 



i-Q 1 
h 



aDlf{t) = lim y (-1 



j=0 



J 



where [a;] means the integer part of x, a,t are the bounds of operation for 
aDUit), g G 3f? and 
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j\{q-j)\ r(j + i)r(g-j + i) 



and 







1. 



We take ideas of a numerical solution method[78, 79] obtained by the 
relationship (11) derived from the G-L definition to resolve system (9) and 
system (10). That is, 

(,_^)Z^Ji/(t)-l^(-l)^ |^^' j/(4_,) 

where is the memory length, tk = kh, h is the time step of calculation and 



J 



are binomial coefficients c^^\ [j = 0, 1, ...,). When for numerical 



computation, the following are used. 



Si) — 1 Ai) 
-0 - -L'S 



1 + ^ M 
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Then in general [78, 79] 
can be expressed as 

fc 

y{tk) = f{y{h),h)h'^-Y,<^fyitk-,) 

j=v 

where v in above is defined as 



V 



IT' ^ IT 



or f = 1 for all k. 

With the ideas above, the method [78, 79] are used here to resolve the 
fractional Lorenz system (9) with a = as following. The general forms of 
chaos systems (12) and (13) are considered respectively as following, where 
i = 1, 77, in equations (12), (13), (14), (15). (/i, /2, fn) in (13) is to be 
reconstructed. 

(12) 

L = (Xi, X2, Xn) 

Xi /i (5^1 5 2^2, •• •, Xn) f-\n\ 

^ ) 

L = (Xi, X2, Xn) 

To have simple forms, we take a = 0. With the ideas above, the systems 
1^12) and (13) are solved as following. 



k 



Xi {tk) = fi {Xi (tfc) , ...,Xi_i {tk) , Xi {tk-l) ,...,Xn {tk-l)) /i"' - ^ cf^'xi (tk-j) 

j=v 

(14) 
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Xi (4) = fi {Xi (tk) , (tk) , Xi (tfc-l) , Xn (tk-l)) h"' - ^ cf'^Xi (tk-j) 

j=v 

(15) 

Then the novel objective function (fitness) equation (16) in this paper 
come into being from equations (14) and (15) as below. 

T-h 2 

f=^||l-l|| (16) 

t=Oh 

The objective function (16) to be optimized can also be any kind of equa- 
tions (3),(4),(5), (6). 

It should be noticed here that the independent variables in function (16) 
are not the parameters but the special variables, for instance, as functions 
(/i,/2, •••,/«)• That is 

(/i,/2,...,/„)* = arg _min_ F (17) 

(/l,/2,...,/n) 

For the reconstructions of three dimensional fractional order Lorenz sys- 
tem (9) and correspond (10), the objective function can be as following. 

/g = argminF (18) 

/a 

Equation (17) and (18) is the crucial turning point that changing from 
the parameters estimation into functions reconstruction, in other words, frac- 
tional order chaos systems' reconstruction. 

Now, with the definitions of equations (16) and (17), it can be concluded 
that the parameters' estimation of fractional order chaos system is a special 
case of fractional order chaos reconstruction. And further, the parameters 
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estimation and reconstruction in Section 2 for the normal chaos system is 
the special cases of fractional order chaos systems' reconstruction (17). 



4. A novel self-growing mechanism of fractional chaos systems' re- 
construction in non-Lyapunov way 

The task of this section is to find a simple but effective approach for 
hyper, proper and improper fractional-order chaos reconstruction. 

4-1. A novel approach for for hyper, proper and improper fractional chaos 
reconstruction 

Now we can propose a self-growing differential equations approach for 
hyper, proper and improper fractional chaos systems. The pseudo-code of 
the proposed reconstruction is given below. 

4-2. Hyper, proper and improper fractional chaos 

To test the Algorithm 1, some different well known and widely used Hy- 
per, proper and improper fractional chaos systems are choose as following. 
Example. 1. Fractional order Arneodo's System (19) [78, 80]. 



^ oDf^z (t) = -Pix (t) - P^y (t) - P^z (t) + (3,x^ (t) . 
when (/3i,/32,/33,/34) = (-5.5,3.5,0.8,-1.0), (gi,g2,g3) = (0.97,0.97,0.96), 



oDrxit)=y{t)- 
< oDry{t) = z{t); 



(19) 



initial point (—0.2,0.5,0.2), Arneodo's System (19) is chaotic. 
Example. 2. Fractional order Buffing's system (20) [78]. 




oDf^x{t) = y{t)- 

QDl'^y it) = x{t)- x^ it) - ay {t) + 5 cos (wt) . 



(20) 



14 



Algorithm 1 Fractional order chaos system reconstruction with genetic op- 
erations 

1: Basic parameters' setting the basic set {x,^, +,— } etc., and the 

input variables {xi,X2, and etc. 

2: Initialize Generate the initial population with tree structures. 
3: while Termination condition is not satisfied do 

4: Evaluation Evaluate the fitness with Eq. (17), where system (13) 
are reconstructed by the individuals. > fractional order differential 

equations 

5: Crossover 

6: Mutation 

7: Reproduction Tournament size. > The Elitism rate is 0.02 — 0.05. 
8: end vi^hile 

9: Output Global optimum (/i, /2, fn) 
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when (a,&,c) = (0.15,0.3,1), (^1,52) = (0.9,1), initial point (0.21,0.31), 
Buffings system (20) is chaotic. 

Example. 3. Fractional order Lotka-Volterra system (21) [78]. 

QDfx (t) = ax (t) - hx {t) y {t) + ex^ (t) - sz (t) (t) ; 
< oDry{t) = -cy{t) + dx{t)y{t)- (21) 
QDfz it) = -pz (t) + sz (t) (t) . 

when (a,6,c,rf,e,p,s) = (1,1,1,1,2,3,2.7), (gi,g2,g3) = (0.95,0.95,0.95), 
initial point (1, 1.4, 1), Lotka-Volterra system (21) is chaotic. 
Example. 4. Fractional order Genesio-Tesi's System (22) [78, 81]. 

oD','xit)=y{t); 
< ,Dry{t) = z{t)- (22) 
^ ^Dfz (t) = -/3ix it) - P^y it) - Psz it) + P,x^ it) . 

when iPi, /32, Ps, P4) = (1.1,1.1,0.45,1.0), (gi,g2,g3) = (1,1,0.95), initial 
point (—0.1,0.5,0.2), Genesio-Tesi's System (22) is chaotic. 
Example. 5. Fractional order Lorenz system (9) with the exact form of 
the differential equation ^Df^z = f^ix,y,z) are unknown, when (a, 6, c) = 
(10,28,8/3), (gi,g2,g3) = (0.993,0.993,0.993), initial point (0.1,0.1,0.1), 
Lorenz system (9) is chaotic. 

Example. 6. Fractional order Lii system (23) [26, 78] . 

' oDrxit) = aiyit)-xit))- 

< oDryit) = -xit)zit) + cyit); (23) 
^ oDrzit) = xit)yit)-bzit). 

when (a, b, c) = (36, 3, 20), (gi, ^2, Qs) = (0.985, 0.99, 0.98), initial point (0.2, 0.5, 0.3), 
Lii system (23) is chaotic. 
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Example. 7. Fractional order Liu System (24) [36, 78, 82]. 

qD^x (t) = -ax (t) - ey"^ (t) ; 
< oDry{t) = by{t)-kx{t)z{t); (24) 

oDl^'z (t) = -cz (t) + mx (t) y (t) . 
when (a, b, c, e, k, m) = (1, 2.5, 5, 1, 4, 4), (gi, q2, gs) = (0.95, 0.95, 0.95), initial 
point (0.2,0,0.5), Liu System (24) is chaotic. 
Example. 8. Fractional order Chen system (25) [37, 78, 83]. 

oD','x{t) = aiy (t) - x (t)) ; 
< oDry (t) = id) x{t)-x it) z (t) + cy {t) ■ (25) 
^ oDrz{t)=x{t)y{t)-bz{t). 
when {a,b,c,d) = (35,3,28,-7), {qi,q2,(l3) = (0.9,0.9,0.9), initial point 
(—9,-5,14), Chen system (25) is chaotic. And when when {a,b,c,d) = 
(35,3,28,-7), (gi,g2,g3) = (1,1.24,1.24), initial point (3.123,1.145,2.453), 
Chen system (25) is an improper chaotic system[20]. 
Example. 9. Fractional order Rossler System (26) [22, 78]. 

' ,Drx{t) = -{y{t) + z{t)); 

< oDry{t)=x{t)+ay{t)- (26) 
^ oDrz{t) = b + z{t){x{t)-c). 
when (a,6,c) = (0.5,0.2,10), (gi,g2,g3) = (0.9, 0.85, 0.95), initial point (0.5,1.5,0.1), 
Rossler System (26) is chaotic. 

Example. 10. Fractional order Chuas oscillator (27) [84]. 

' ^Dfx (t) = a{y{t)-x (t) + (x (t) -W{w)x (t)) ; 

oDry{t) = xit)-yit) + z{t); ^^^^ 

^Dfz{t) = -l3y{t)-iz{t)- 
^ oDpw{t)=x{t); 
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where 



W{w) 



1 a : 


\w\ 


I- 


\w\ 



(28) 



when {a, (3, 7, C, a, b) = (10, 13, 0.1, 1.5, 0.3, 0.8), (gi, ga, gs, ^4) = (0.97, 0.97, 0.97, 0.97), 
initial point (0.8,0.05,0.007,0.6), Chua's oscillator (27) is chaotic. 
Example. 11. Hyper fractional order Lorenz System (28) [85] . 

oDl'x{t) = a{y{t)-x{t))+w{t); 
oDfy (t) = cx (t) -x{t)z (t) -yit); 
oDrz{t)=x{t)y{t)-bz (t); 
oDTw it) = -y it) z it) + -iw it) ■ 

when {a,b,c,d) = (10,8/3,28,-1), (gi, ^2, gs, ^4) = (0.96,0.96,0.96,0.96), 
initial point (0.5,0.6,1,2), Hyper fractional order Lorenz System (28) is 
chaotic. 

Example. 12. Hyper fractional order Lii System (29) [86]. 

oDl'x{t) = a{y {t) - x {t)) + w {t) ■ 
oDry{t) = -x {t)z{t) + cy (t) ; 
oDrz{t)=x{t)y{t)-bz (t); 
qD^w it) =x{t)z it) + dw it) ; 

when (a, 6, c, d) = (36, 3, 20, 1.3), (gi, ga, gs, q^) = (0.98, 0.980.98, 0.98), initial 
point (1, 1, 1, 1), Hyper fractional order Lii System (29) is chaotic. 
Example. 13. Hyper fractional order Liu System (30) [87]. 

oA^'x (t) = -ax (t) + by (t) z (t) + z (t) ; 
oDry{t) = 2.5y (t) - x (t) z (t) ; 
oDf'z (t) = x{t)y (t) - cz (t) - 2w (t) ; 
oDpw (t) = -d-x{t). 



(29) 



(30) 
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when {a,b,c,d) = (10,1,4,0.25), (gi, g2, gs, ^4) = (0.9,0.9,0.9,0.9), initial 
point (2.4,2.2,0.8,0), Hyper fractional order Liu System (30) is chaotic. 
Example. 14. Hyper fractional order Chen System (31) [88]. 

oDl'x{t) = -a{y{t)-x{t)) + w{t); 
oDTy{t) = dx{t)-x{t)z{t) + cy{t)- ^^^^ 
,DTz{t) = x{t)y{t)-hz{t)- 
^ oDl'w (t) =y{t)z (t) + rw (t) . 

when {a,b,c,d) = (35,3,12,7,0.5), (gi, ^2, gs, ^4) = (0.96,0.96,0.96,0.96), 
initial point (0.5, 0.6, 1, 2), Hyper fractional order Chen System (31) is chaotic. 
Example. 15. Hyper fractional order Rossler System (32) [22]. 

' oDrx{t) = -{y{t) + z{t)); 
oDry{t)=x{t) + ay{t) + w{t)- 

(32) 

oDTz{t)=x{t)z{t) + h- 
oD^*w (t) = -cz (t) + dw (t) . 

when {a,b,c,d) = (0.32,3,0.50.05), (gi, g2, gs, g4) = (0.95, 0.950.95, 0.95), ini- 
tial point (—15.5,9.3, —4, 18.6), Hyper fractional order Rossler System (32) 
is chaotic. 

4-3. Simulations 

There is one important thing we have to deal with is that the normal ODE 
methods are not valid for many individuals in evolution of CP. Even they 
are efficient in some cases, however the time consuming is too much. Here 
we assign all the values of NAN and INF as for fractional chaos systems' 
resolving, then the process of getting results by resolving the fractional order 
system can be more simpler. 
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For systems to be reconstructed, the parameters of the proposed method 
are set as following. Population size 100, max generations 200, mutation 
rate 0.05, crossover rate is 0.85, directed copy rate is 0.02 (the Elitism rate), 
the tournament selection is implemented (the size of tournament is 3), the 
maximum depth of tree is 6, the maximum depth of subtrees created by 
mutation is 3, the maximum number of nodes per tree is infinite, and the 
input variables for 3 and 4 dimensional fractional order chaos systems are 
taken and x,y,z,w respectively. Table 1 give the detail setting for 

each system. 

Table 2 shows the simulation results of above fractional order chaotic 
systems. 

The following figures give a illustration how the self growing evolution 
process works by genetic operations of Algorithm 1. In which. Figures 1,2 
,3, 4, 5, 6,7 ,8, 9, 10 show the simulation evolution results of above fractional 
order chaotic systems. Figures 11,12 show some of the correspond objective 
function's evolution process of above fractional order chaotic systems. 

From the simulations results of reconstruction above fractional order 
chaos system, it can be concluded that the proposed method is efficient. 
And from above figures, it can be concluded that the estimated systems are 
self growing under the genetic operations of the proposed methods. 

5. Conclusions 

The put method consists of numerical optimization problem with un- 
known fractional order differential equations to identify the chaotic systems. 
Simulation results demonstrate the effectiveness and efficiency of the pro- 
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Table 1: Detail parameters stettiiig for different systems 



F-0 system 


Unknown 


Basic set 




Constant 


Probabihty 


Step 


No. 


Arneodo 






-} 


[-5,5] 


0.1 


0.005 


200 


Duffing 


D^x 


{x,+,-} 








0.0005 


500 


Lotka-Volterra 


D^x 


{x,+,-} 








0.01 


100 


GenTesi 


Dfx.DTy 


{x,-,+,- 


-} 


[-5,5] 


0.1 


0.005 


200 


Improper Chen 


Dfx 


{x,+,-} 




[-3,3] 


0.3 


0.005 


200 


Lorenz 


Dfx 


{x,+,-} 








0.01 


50 


Lorenz 


Dfz 

t 


{x,+,-| 




[-3,3] 


0.5 


0.01 


50 


Lii 


DTz 


{x,+,-} 




[-5,5] 


0.2 


0.01 


100 


Liu 


Dfx 


{x,+,-} 








0.01 


100 


Chen 


Dfz 


{x,+,-} 




[-5,5] 


0.2 


0.01 


100 


Rossler 


DTz 


{x,+,-} 








0.01 


100 


ChuaM 


D?y 


{x,+,-} 








0.01 


100 


Hyper Lorenz 


Dfw 


{x,+,-} 








0.01 
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Hyper Lii 


DTz 


{x,+,-} 








0.005 
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Hyper Liu 


Dfw 


{x,+,-} 




[-1,1] 


0.5 


0.005 


100 


Hyper Chen 


Dfz 


{x,+,-} 








0.005 


200 


Hyper Rossler 


DTz 


{x,+,-} 




[-4,4] 


0.5 


0.005 


200 



of samples 
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Tabic 2: Simulation results for different fractional order chaos systems 



F-0 system Unknown StD Mean'* Min Max Success rate^ 
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yyzz.o ( zo 


00/0 


T 


J->t ^ 




Z.iDDDe-4 


o.y4oie-o 


n 
U 


i.Dyzoe-o 


Q1 (J/ 

yi/o 


Chen 






4oo.00oo 


1 /ic nom 
i4D.yzuz 


n 
U 


ofinn /lOQO 
Z0UU.4Z0Z 


< U/o 


Rossler 
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17742.4778 


79% 


Hyper Lii 


Dfz 




1263.4313 
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95% 
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56% 
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Dfz 




4754.4523 
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0.8458 


0.1373 


4.8299e-14 


6.0346 


62% ' 



Std and Mean are calculated for all the solutions. 



Success means the the solution is less than le — 20 in 100 independent simulations. 
Success means the the solution is less than le — 1 in 100 independent simulations. 
Success means the the solution is less than le — 8 in 100 independent simulations. 
Success means the the solution is less than le — 4 in 100 independent simulations. 




Figure 1: The fractional Chen system with best Df^ z 
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Generation 9 

^ {b)Dfx{t) =1/ ~yz~x^ ~x + y^ 

Generation 1 

{a)Dfx{t) = z-x-{x-y + xz{y + z)){2x + 3y + z) -xz-{x + y)(y - z] - 2z' 

2 . 




Generation 10 Generation 15 

(c)D'j' X it) = -r' -x + 'if -yz {d)Df x{t) = if - 



Figure 2: The fractional Liu system with best Df^x 
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Generation 1 

(a)r)fx(t) =2y + 2z 



Generation 16 



60^ 

40, 

20, 

0, 
50 




-50 



20 



40 



-2% 

Generation 36 

{c)Dfx(t) = ly-2x-2z 




Generation 46 

(d)Dfx(t) = %y-ix-2z 



40, 
20, 




-50 _20 -10 ,. 10 20 30 



Generation 66 

(e)Dfx(t)=9y-7x-z 




10 20 30 



Generation 81 
U)Dfx(t)=10y-Wx 



Figure 3: The fractional Lorenz system with best D'^^x 
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Figure 4: The fractional Lii system with best Dpz 
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Figure 5: The fractional Rossler system with best Df^ z 
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-10 a; 10 20 , 
Generation 1 
(a)Dfz(t)=xy 




-20 X 

Generation 9 
(b)Dfz{t)=Zy- 



-xy 




-4 -nZO 20 w 
Generation 16 

{c)Dfz(t) = 'Sy-x-2z + xy 




-40 ^0 20 w 
Generation 32 

{d)Dfz(t) =iy-2x-2z + xy 
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30 
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„ 20 
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50 



-j2 2 lii 
Generation 80 

(e)Dfz (t) = :3y-'ix-2z + xy 




-20 



Q. 2 4 u. 
Generation 103 
(J)Dfz{t)=xy-Zz 



Figure 6: The hyper fractional Chen system with best Df'z 
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Generation 40 x 

(c)Dfw (t) = (-0.022042^ - 0.2435)x Generation 50 

{d)Dfw{t) = (-0.014442^ - 0.2455)x 




Generation 110 Generation 140 

{e)Efi'w (i) = (-0.01235Z - 0.2435)i {DD'l'w (i) = -0.252: 



Figure 7: The hyper fractional Liu system with best Dl^w 
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Figure 8: The hyper fractional Lorenz system with best D^^w 
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Figure 9: The hyper fractional Lii system with best D^^ z 
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Generation 4 Generation 15 

{a)Dfz [t) =xz- 3.089Z + 0.09988 (b)Dfz (i) = xz + 2.583 




Figure 10: The hyper fractional Rossler system with best Df^ z 



32 




Figure 11: Evolution process of objective function F 
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Figure 12: Evolution process of objective function F 
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posed methods with the unified new mathematical modeL This is a novel 
Non-Lyaponov way for chaos reconstruction. 

We still have to improve the searching performance for genetic opera- 
tions of proposed method, for the reason that the more efficient the genetic 
operations are, the better the new method is in some degree. Moreover, it 
also illustrated the simplicity and easy implementation of the other kind of 
intelligent algorithms to replace genetic operations of Algorithms 1 for appli- 
cations, as in genetic operations from differential evolution algorithms added 
to tree structures in reference [89]. 

The performance of the proposed method is sensitive to the initial point 
for each fractional order chaos system, sample period, number of points, and 
intervals. Actually, these also lead to the candidate system divergent. And 
they are not predefined randomly. A good combination of these is not easy 
to get. Some mathematical formula to get a good combination not by so 
many simulations will be introduced in the future studies. 

It should be noticed too many points for evaluating chaos system the 
individual represents are not worth. Because the most time consumption 
parts in the whole proposed method are to resolve the candidate systems. 
Some of these system are easy to solve. However when it comes with the 
bad individuals, the methods to resolve the fractional order chaos systems in 
Section 3 might not converge. Then the whole proposed method Algorithms 
1 might get into endless loops. To avoid the endless loops, we introduce a 
forced strategy to assign all the NAN and infinite numbers in the output as 
zero. Because the objective function (17) to be optimized is bigger than 0, 
so this forced strategy for assignment is reasonable. To achieve a fine bal- 
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ance between the performance of the proposed methods and having enough 
sample data for credibihty, we take the number of the points as 100 — 200, 
according the existing simulations [40, 41, 43-45, 47, 48, 50, 51, 90-95]. And 
the simulations in section 4 results show it is effective too. 

Here we have to say the success rate has something to do with the prob- 
lems to be reconstructed themselves. No methods are valid for all the prob- 
lems. According the success rate of Table 2, there are different results for 
different fractional chaos systems, even the same methods are used. In which, 
the worst one is for fraction order Lorenz system and hyper fractional order 
Rossler and Liu systems. The reason lies in the systems themselves. For Df^z 
of the fraction order Lorenz system, the best one is D^^ z = x ■ y — b ■ z with 
b = 8/3. However, it is not easy to achieve 6 = 8/3. As we know from the 
existing simulations [40, 41, 43-45, 47, 48, 50, 51, 90-95], even the estimated 
b is very near to real value 

b = 8/3 

, the estimated system are not the same as the real Lorenz system. Thus for 
the Algorithms 1 here, it has two tasks, not only finding the best form for 
fractional differential equations Df^z = f^, but also finding the best value for 
b. For the hyper fractional order Rossler system, Df^z = xz + 3.0 is not easy 
to be achieved either. Here the key point is the constant 3 G D^^z = xz + 3.0. 
And the same to the hyper fractional order Liu system where 0.25 G Dl'^w = 
0.25x 

In conclusion, it has to be stated that self growing fractional order dif- 
ferential equations in chaos systems' reconstruction is a promising direction. 
And the other EAs can also be introduced inside of the proposed methods, 
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by doing their own special genetic operations to the individuals in GP with 
the objective function (17) . 
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